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Abstract 

In this paper we apply a well-tested approximation of electron Coulomb 
distortion effects to the exclusive reaction (e, e'p) in the quasielastic region. 
We compare the approximate treatment of Coulomb distortion effects to the 
exact Distorted Wave Born Approximation evaluated by means of partial 
wave analysis to gauge the quality of our approximate treatment. We show 
that the approximate Mpller potential has a plane-wave-like structure and 
hence permits the separation of the cross section into five terms which depend 
on bi-linear products of transforms of the transition four current elements. 
These transforms reduce to Fourier transforms when Coulomb distortion is 
not present, but become modified with the inclusion of Coulomb distortion. 
We investigate the application of the approximate formalism to a model of 
208 Pb(e, e'p ) using Dirac-Hartree single particle wavefunctions for the ground 
state and relativistic optical model wavefunctions for the continuum proton. 
We show that it is still possible to extract, albeit with some approximation, 
the various structure functions from the experimentally measured data even 

for heavy nuclei. 
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I. INTRODUCTION 


Medium and high energy electron scattering has long been acknowledged as a useful tool 
in the investigation of nuclear structure and nuclear properties, especially in the quasielastic 
region. In the plane wave Born approximation (PWBA), where electrons are described by 
Dirac plane waves, the cross section for the exclusive reaction (e, e'p) on nuclei can be written 
simply as 
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where q 2 = oj 2 — q 2 is the four-momentum transfer, % is the Mott cross section given 
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f, and Rp, Rp, Rtt, Rlt, and Rlt 1 are the longitudinal, transverse, 


transverse-transverse, longitudinal-transverse, and polarized longitudinal-transverse struc¬ 
ture functions which depend only on the momentum transfer q and the energy transfer uj. 
The functions v only depend on electron kinematics and are given by 
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Choosing the momentum transfer q to define the z axis and using the azimuthal sym¬ 
metry of the spatial part of the Mpllcr potential permits the extraction of the explicit 
dependence on the azimuthal angle 0p of the outgoing proton as measured with respect to 
the electron scattering plane. More specifically, we define y = p, x pj. These structure 
functions are defined as: 
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where the nuclear tensors are given in terms of a sum over s p and /i p the final and 
initial spin projections of the nucleon, 

nv =e n ; n - < 4 ) 

Spflb 

and we have suppressed the spin labels for clarity. The quantities iV M are the Fourier 
transform of the nucleon current density J /i (r) given by 

AT = J (5) 

Current conservation and gauge invariance can be used to eliminate the ^-components so 
that only A^, N x , and N y need to be calculated. 

By keeping the momentum and energy transfer fixed while varying the electron energy 
E and scattering angle 9 e , or varying the azimuthal angle of the outgoing proton and the 
helicity of the incident electron, it is possible to extract from experiment the five structure 
functions as a function of momentum and energy transfer. However, when the electron 
wavefunctions are not Dirac plane waves, but rather are distorted by the static Coulomb 
held of the target nucleus, such a simple formulation as given in Eq. (1) is no longer possible. 
In the full Distorted Wave Born Approximation (DWBA) calculation, it is not possible to 
express the cross section as a sum of bi-linear products of transforms of transition current 
matrix elements which only depend on the electron kinematics and the outgoing proton’s 
azimuthal angle. The key point is that the momentum transfer does not enter the analysis 
in a natural way. Of course, one could pretend that the plane-wave result is still valid and 


3 



extract the various ”structure functions’. However, there is no way in the DWBA approach 
to investigate these terms separately. 

Using partial wave analysis, the Ohio University group [1]- [4] has treated the Coulomb 
distortion arising from the static Coulomb held of the nucleus exactly. Coupling the dis¬ 
torted waves with the one hard photon exchange approximation (Distorted Wave Born 
Approximation-DWBA) has allowed this analysis to be compared to data from a range of 
nuclei. The nuclear model used includes the following ingredients: 1) Relativistic Hartree 
single particle wavefunctions for the bound orbitals [5,6], 2) Relativistic optical model for 
the continuum proton [7], and 3) the free relativistic current operator for the proton with 
the standard form factors. This simple relativistic “one-body” model along with the ex¬ 
act treatment of Coulomb distortion is in excellent agreement with all the data involving 
knock-out from surface orbitals of nuclei that have been analysed. This includes nuclei as 
light as 16 0 and as heavy as 20S Pb. However, for cases where the outgoing particle is in 
the continuum where all multipoles can contribute, the DWBA analysis requires extensive 
computer codes which require more and more time and precision as the energy increases. 
Furthermore, as noted above in the exact DWBA analysis, the cross section cannot be writ¬ 
ten as the sum of five terms which are bi-linear products of transforms of the transition 
current matrix elements. There two reasons have led researchers to seek an approximate 
treatment of Coulomb distortion that would permit the extraction of “structure” functions 
from experiment under conditions where the effects of Coulomb distortion in the different 
terms can be investigated, and could be easily extended to higher energies such as will be 
available at the Thomas Jefferson National Accelerator Facility. 

In a previous paper [8] we investigated a rather extreme approximation of Coulomb 
distortion effects for the inclusive reaction (e, e') in the quasielastic region and using some 
ad — hoc assumptions obtained excellent agreement with the exact DWBA calculations. 
In this paper, we apply a more exact approximation to the exclusive reaction (e, e'p) and 
investigate its validity. One of our advantages as compared to previous investigators is that 
we have the full DWBA calculation to use as a standard in assessing the accuracy of our 
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approximation. In Section II we will give the approximate electron wavefunctions for the 
incoming and outgoing electron waves and derive the approximate Mpllcr potential from 
these wavefunctions. In Section III we use the approximate Mpller potential to calculate the 
differential cross section for (e, e'p) and define the approximate Coulomb deformed structure 
functions. In Section IV we compare approximate and DWBA (e, e'p) cross sections for so- 
called parallel and perpendicular kinematics and investigate the extraction from the cross 
sections of the so-called fourth (Rlt) and fifth ( Rlt') structure functions. In Section V we 
give our conclusions and discuss prospects for the future. 

II. APPROXIMATE WAVEFUNCTIONS AND THE M0LLER POTENTIALS 

Following the work of Lenz and Rosenfelder [10], we propose the following plane-wave- 
like electron wavefunction [8,9] which contains the effects of the static Coulomb distortion 
of the target nucleus in an approximate way: 

T (±) (r) = ^i e ±w(j2) e iA e <p ' (r) - r « p . (6) 

P 

The (±) sign denotes incoming and outgoing boundary conditions for electrons of momentum 
p, the phase factor <5(J 2 ) is a function of the square of the angular momentum operator J, 
and the local effective momentum p'(r) is given in terms of the Coulomb potential of the 
target nucleus by 

p'(r) — (p—~ [ V(r)dr) p. (7) 

r Jo 

We refer to this r-dependent momentum as the Local Effective Momentum Approxi- 
mation(LEMA). Some small higher order corrections have been incorporated into the 
ad — hoc term A = a(p , (r)-r)(J 2 + |) which involves the factor a which is parametrized 
by a = — aZQhJ) 2 where Z is the charge of the target nucleus and the number 16 is given in 
MeV/c and was determined by comparison with the exact result. Inclusion of this term with 
LEMA is referred to as LEMA + A. The electron mass has been neglected in comparison to 
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the electron momentum, so this approximation if not valid at extreme forward and backward 
electron scattering angles. 

Previous workers [12,13] replaced the r-dependent momentum p'{r) in Eq. (7) by the value 
Pema = p~ V(0). This approximation is known as the Effective Momentum Approximation 
(EMA). Unfortunately, as we have known, the EMA describes the Coulomb effects on the 
wavefunction rather poorly [8,9]. Previous workers [12,13] also approximated the Coulomb 
phases by a constant plus a linear term in the operator J 2 . While this approximation works 
well for low partial waves, it does not describe partial waves with angular momenta equal 
to or greater than pR where R is the nuclear radius. However, it is these partial waves that 
dominate inelastic electron scattering from the nucleus. We refer to the EMA plus linear 
fit to the phases as EMA-k 2 . We avoid this problem by fitting the exact partial wave phase 
shifts S K , where k is the Dirac quantum number, to an expansion in powers of k 2 . Retaining 
terms of second order in n 2 for kappa values up to approximately 3 (pR), we fit the exact 
phases with the following equation, 

S K = b 0 + b 2 K 2 + (8) 

Note that the eigenvalues of J 2 are j (j + 1) which equals k 2 — 

We also investigated the expansion of the phase exponential in Eq. (6) into a power 
series carried out by previous workers [12,13] and concluded that an accurate description 
requires many terms. We chose instead to neglect the electron spin dependence of the phases 
and replace J 2 + 1 by the angular momentum squared L 2 in the exponential phase term for 
both incoming and outgoing waves. Further, we replace L 2 by its classical value (rxp') 2 . 
With these two approximations, the approximate Coulomb wavefunction is given by 

V|/(±)( r ) — Pjy}_ c ±ib 0 c ±ib 2 {rxp(r )) 2 

P 

xe ±ib 4 (rxp(r)) 4 e ia(p , (r)-f)(rxp(r)) 2 e jp'(r)-r M ^^ (g) 

The merits of this approach are that our approximate distorted Coulomb wavefunctions have 
an analytic plane-wave-like form and the Coulomb distortion modifications do not depend 
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on the electron spin, but the wavefunctions do have (rxp'(r)) 2 terms in the exponential 
which carry information about the phase shifts. We will refer to this wavefunction as the 
approximate analytic Coulomb distorted wavefunction. Based on our investigations, it is a 
good representation of the exact partial wave solutions for radial coordinates out to about 
three or four nuclear radii. 

Using a technique introduced by Knoll [11] which approximates the potential in terms of 
the source current we obtain the following approximate Mpller-type potential corresponding 


to our approximate analytic Coulomb distorted wavefunction, 
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The Knoll approach is discussed in detail in previous work [8,9] and is a good approximation 
for momentum transfers greater than about 350 MeV/c. In arriving at Eq. (10), we neglected 
spatial derivatives of the phase factors in the wavefunction and the radial derivative of the 
local effective momenta p'(r). 

The comparison of the approximate four potentials with the DWBA four potential re¬ 
quires a partial wave expansion since that is the only way the DWBA results can be obtained. 
However, the partial wave expansion for these potentials is hardly possible because of the 
(rxp') 2 terms in the exponential. Traini et al., [13] expanded the exponential in a power 
series up to the second order so as to use partial waves, but this series converges very slowly 
because the value of the phase, 6(rxp) 2 , is greater than one for regions of space near the 
nuclear surface. Of course, the origin of our approximations was in a partial wave formalism, 
so we can go back to that formalism and replace the exact phases and radial wavefunctions 
by the approximate phases and wavefunctions to obtain a measure of the quality of our 
approximations. Note, however, that obtaining the plane-wave-likc form required some ad¬ 
ditional approximations including assuming that the asymptotic momentum transfer q and 
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the phase factors in the approximate four potential do not explicitly depend on the incoming 
and outgoing electron spins. 

With this caveat, we can compare various approximate potentials that have been widely 
used [12,13] and our approximate potential for selected multipoles. Figs. (1) and (2) show 
the comparison of two scalar potentials with the DWBA scalar potential by using the partial 
waves for different multipole L values. The initial spin and the final spin of the electron 
are s l = Sf = 1/2. The calculations have been done using a Fermi charge distribution of 
radius R = 6.65 fm, with total charge Z = 82 and the angular momentum L = 5 and 
15. The initial electron energy is Ei = 400 MeV, the final electron energy is Ef = 300 
MeV, and the momentum transfer is q = 350.5 MeV/c. The solid line is the result for the 
DWBA scalar potential, the dotted line is for the EM A — E 2 potential, and the dashed 
line is for our approximate potential of Eq. (10). The approximation for EM A — k 2 has a 
different magnitude and is also out of phase, but that for our approximate potential almost 
has the same magnitude and is in phase. The discrepancy at large radii is due to lack 
of complete convergence of the n series in the DWBA calculation and does not play an 
important role in calculating the cross section since the bound nucleon wavefunction drops 
very rapidly with radial distance. We conclude that our approximate potential is in quite 
good agreement with the exact potential calculated with partial waves (DWBA) for radial 
coordinates less than three to four nuclear radii, and hence is good enough to replace the 
full DWBA calculation. We have confirmed with the aid of a simple model [9] that our 
approximate potential reproduces the cross section calculated with the full partial wave 
result from DWBA quite well. We will show more realistic comparisons in the following 
sections that do not utilize a multipole decomposition and therefore is a more direct test of 
the approximate Coulomb distorted potential given in Eq. (10). 



III. APPROXIMATE COULOMB DISTORTED CROSS SECTION 


Using the approximate Mqllcr potential given in Eq. (10) it is straightforward to derive 
the cross sections for (e, e'p) reactions from nuclei since apart from the modified spatial 
dependence the approximate potential has the same Dirac structure as the plane wave Mpller 
potential. The result is 
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The electron structure functions are unchanged from Eq. (2), but the nuclear structure 
functions, designated with a prime superscript, contain Coulomb distortion effects. They 
are defined as in Eqs. (3) and (4), except that the “Fourier” transforms given in Eq. (5) 
become 
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where the phase shift 5 and the ad-hoc additional phase A are functions of r given by 


5 = b 0 + & 2 (rxp(r)) 2 + 6 4 (rxp'(r)) 4 

A = a(r-p(r))(rxp'(r)) 2 . (13) 


The result in Eq. (11) is our primary finding. Our approximate treatment of Coulomb 
distortion leads to a “plane-wave-like” form for the cross section and thereby opens up 
the possibility of investigating the various ’’structure functions” independently. Of course, 
these more generalized “structure functions” do contain some dependence on the electron 
kinematics, but with the use of theoretical nuclear models, this modification of the structure 
functions can be investigated. 

The charge transform in Eq. (11) differs from the transverse current transforms since 
the continuity equation was used to eliminate the z-component of the current. Note that 
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unlike the case for electron plane waves, the various current transforms are not azimuthally 
symmetric about the momentum transfer direction q, and therefore contain dependencies on 
the outgoing nucleon azimuthal angle (pp over and beyond the explicit dependencies shown 
in Eq. (11). However, some symmetry remains since both (rxp;'(r)) 2 and (rxpf'(r)) 2 are 
invariant under the transformation <p —► —0 which results in the nuclear structure being 
invariant under <pp —> —0p. The consequences of this additional dependence on (pp will be 
discussed further in the next section. 


IV. RESULTS 


A. Cross section in parallel and perpendicular kinematics 


In our analysis we are looking at one particular shell, and trying to find the reduced 
cross section p m , which for proton waves in the final state is related to the probability that 
a bound proton from a given shell with the missing momentum p m can be knocked out of 
the nucleus with asymptotic momentum P. The reduced cross section as a function of p m 
is commonly defined by 


Pm(p m ) PE P a eP dE f dfl f dfl P ’ ^ 

where the missing momentum can be determined by the kinematics p m = P — q. The 
off-shell electron-proton cross section a eP is not uniquely defined, but in all cases we use the 
form (Tgp given by deForest [14]. 

There are two kinematical situations commonly used in (e, e'p) experiments. They are 
parallel kinematics where the outgoing proton momentum P is along the momentum transfer 
q and perpendicular kinematics where the magnitude of P is fixed, but the detected proton 
makes an angle 6p q with respect to q. In perpendicular kinematics, the magnitude of P is 
usually equal to the magnitude of q. All calculations will be carried out in the laboratory 
frame (target fixed frame). In the parallel case, the three interference terms in Eq. (14) dis¬ 
appear, while in the perpendicular case, all terms remain except the fifth structure function 
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which sums to zero for unpolarized incoming electrons. Our approximate calculations for the 
(e, e'p) reaction include the approximate phase factors and the correction term A by keeping 
the exponential form and the transition matrix element is evaluated by three dimensional 
integration since a multipole expansion is no longer practical. We compare our results to the 
full DWBA results [1,2] and the experimental data from NIKHEF [15,16] in Amsterdam. 
The electron incoming energy is given by E t = 412 MeV and the ejected proton kinetic 
energy Tp = 100 MeV. All calculations include the proton final state interaction by using a 
relativistic optical potential obtained from fitting to elastic proton scattering data [7]. 

In Figs. (3) and (4), we show two results corresponding to knocking out a proton from a 
3 si/ 2 shell and a 2d 3 / 2 shell in 208 Pb for the case of parallel kinematics. For this kinematics, 
the proton momentum P is parallel to the asymptotic momentum transfer q which defines 
the z axis and the missing momentum p. m is also along the q direction. The dotted line is the 
PWBA result obtained by using a multipole expansion and doing the one dimensional inte¬ 
gration over r in the normal way, while the dash-dotted line uses the approximate potential 
with Z=1 evaluated by three dimensional numerical integration. They are in excellent agree¬ 
ment (to better than 1%) as they should be since the approximate calculations approaches 
the plane wave result as Z— >0. The solid line is our approximate Coulomb distorted result 
obtained by numerical integration while the dashed line is the full DWBA results obtained 
by using partial waves and multipole analysis. The dash-three-dotted line is the EM A — E 2 
result also obtained by using three dimensional integration. The diamonds are data from 
NIKHEF. Note that the primary effect of Coulomb distortion is to shift the reduced cross 
section by about 20 MeV/c in missing momentum to the right as compared to the plane 
wave results. The approximate DWBA results reproduce the full DWBA results well around 
the first peak where the difference is less than 2%, but deviate somewhat for large miss¬ 
ing momentum where the reduced cross section is smaller. The approximate DWBA result 
breaks down rapidly beyond missing momentum p m = 100 MeV/c on the right side, but 
since q = P — p m , positive p m corresponds to small q and we expect our approximation to 
become worse for q less than about 350 MeV/c as discussed in Section 2. 
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The EM A — k 2 result is lower by about 30% around the first peak than the full DWBA 
result. Note further that the electron distortion affects the positive missing momentum p m 
and the negative p m differently. The negative p m region shows a large Coulomb distortion 
effect. Figs. (5) and (6) show the reduced cross sections from 3si/ 2 and 2d 3 / 2 for 20S Pb 
for perpendicular kinematics. We choose P = q which marks the top of the quasielastic 
peak for s-states. The electron angle corresponding to this momentum transfer is 6 e = 74° 
for Ei = 412 Mev and the ejected proton energy Tp = 100 MeV. The dotted line is the 
PWBA result, the solid line is our approximate result, the dash-three-dotted line is the 
EM A — k 2 result, and the dashed line is the full DWBA result. Since the momentum 
transfer is large, q = 444 MeV/c, the approximate results are in much better agreement 
with the full DWBA results than those of the parallel kinematics case discussed above. The 
difference is less than 2% around the first peak as in the parallel case, and for both side 
regions the deviation is around 5%. The positions of the maxima and minima are in the 
right places. The discrepancy between the DWBA result and the EM A — k 2 result is greater 
than 30% around the first peak. These results confirm previous observations that Coulomb 
distortion has smaller effects in perpendicular kinematics than in parallel kinematics. 

From these calculations of two different kinematic cases for the (e, e'p) reaction, the 
approximate DWBA results reproduce the full DWBA and the experimental results quite 
well, especially around the first peak. The effect of the Coulomb distortion on the cross 
section for a knocked-out proton from the 3si/ 2 shell for 208 Pb is almost 30% but that for 
knocked-out from the 2d 3 / 2 shell is only 10% as compared to the PWBA calculation. 

In the past, (e, e'p) experiments in parallel kinematics have been measured for the re¬ 
duced cross section in the missing momentum range —50<p m <300 MeV/c at NIKHEF. 
Recently, the range of the missing momentum has been extended to 300<p m <500 MeV/c 
by (e, e'p) measurements for perpendicular kinematics [17]. The new reduced cross section 
was measured at momentum transfer q = 221 MeV/c, energy transfer u = 110 MeV, the 
kinetic energy of the detected proton T p = 100 MeV, and incident electron energy E t = 487 
MeV as shown in Fig. (7). The dotted line is the PWBA result, the solid line is our ap- 
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proximate DWBA result, and the diamonds are the experimental data from NIKHEF [17]. 
Our result reproduces the measured reduced cross section very well although the momentum 
transfer q is small. A similar conclusion has beed drawn [19] by the Madrid group using 
their DWBA calculation to analyse this same data. The interesting physics point is that 
our “single-particle” relativistic model reproduces the experimental data at large missing 
momentum quite well and should be compared to an analysis of this same reaction with a 
non-relativistic approach which uses a non-relavistic current operator and finds it necessary 
to introduce many different two-body currents to even come close to the data [18]. It should 
be noted that our calculation only contains one free parameter, the spectroscopic factor 
which is an overall scale factor of 0.71 which had already been determined [1]- [4] by the 
low missing momentum data. 

B. Interference Structure Functions 

In previous work [4], the effect of Coulomb distortion on the magnitude of the fourth 
structure function was more than 15% in 16 0, and more than a factor of 2 in 208 Pb. The 
magnitude effect on the fifth structure function depends on the out-of-plane angle of the 
knocked out proton used in the extraction and was more than 15% for a small angle (e.g. 
10°) data in 16 0. Since the structure functions appear in the cross section with different 
electron kinematic factors, one can study them independently, but we will show below that 
the PWBA formalism is no longer valid in the presence of the static Coulomb held of 
the nucleus. Even though the separation for the full DWBA calculation with a partial 
wave expansion is not valid in the presence of the Coulomb distortion, it is possible to 
calculate the fourth and the fifth structure functions which embody left-right and up-down 
asymmetries of the cross section measured with respect to the momentum transfer direction. 
We call a quantity so determined the apparent structure function, and note that it would 
correspond to a structure function extracted from experiment. In our model, we can also 
directly calculate the ”structure functions” as given by Eq. (3) when the distorted ’’Fourier” 
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transforms of Eq. (12) are used. One question is to what extent these two results agree with 
each other. 

From Eq. (1) one can see that the fourth structure function Rlt could be obtained 
experimentally by subtracting the cross sections with (j)p = 0 and (f>p = n while keeping 
other electron and proton kinematic variables fixed. The fourth apparent structure function 
determined by the left-right asymmetry with respect to the momentum transfer direction is 
given by 


TDCL _ 

SX L T — 


a L _ fjR 


2 Kv 


LT 


(15) 


where L (left) indicates (f)p = 0 and R (right) indicates <f)p = n. The constant K = 
and the electron structure functions v are defined in Eq. (2). The superscript a means the 
apparent structure function including the electron Coulomb distortion, and corresponds to 
what one would extract from experiment. If the incoming electron beam is polarized (h = 1), 
one can extract the apparent fifth structure function by the up-down asymmetry of the cross 
section given by 
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2ATy T /sin (f) P 


(16) 


where U (up) means 0 < (pp < n (above the plane) and D (down) means —<pp (below the 
plane). 

We extract the fourth structure functions for the 3 Si/ 2 orbit of 208 Pb with incident elec¬ 
tron energy E t = 500 MeV as shown Fig. (8). The dotted line and the dashed line are the 
fourth structure functions calculated directly using Eq. (12) for the PWBA and the approxi¬ 
mate DWBA result, while the dash-dotted line, the solid line, and the dash-three-dotted line 
are the apparent fourth structure functions from Eq. (15) for the PWBA, the approximate 
DWBA and the full DWBA calculation, respectively. We compared our approximate DWBA 
calculation of the cross section to the full DWBA and found the difference to be less than 
2% around the first peak, and around 5% at the second peak. When Coulomb distortion is 
included the directly calculated fourth structure function differs from the apparent fourth 
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structure function by a factor of 3 at the peaks. Of course for the PWBA calculation the ap¬ 
parent fourth structure function agrees exactly with the directly calculated fourth structure 
function. Clearly, the standard separation formalism is no longer valid in the presence of 
the electron Coulomb distortion for the fourth structure function Rlt- Furthermore, while 
the effect of the electron Coulomb distortion is on the order of 30% for the cross section, it 
changes the apparent fourth structure function by more than a factor of 2. 

In Fig. (8) the discrepancy between the directly calculated fourth structure function 
R' lt and the apparent fourth structure function for the approximate DWBA suggests that 
the structure functions depend on the azimuthal angle of the ejected proton as expected. 
In order to reduce this dependence, we investigated changing the definition of the z axis, 
normally defined by asymptotic momentum transfer q, in order to bring the apparent and 
direct structure functions into closer agreement. We considered two choices, one where z is 
taken to lie along q \R) = p [(R) — p 'f(R) where p'(R) = p — V(R) and V(R) is the value of 
the Coulomb potential at the nuclear surface, and the second along q'(0) = p'(0) — p^(0) 
where p'(0) = p — V(0) and V(0) is the Coulomb potential at the origin. The second 
case corresponds to the EMA approximation. We carried out our approximate DWBA 
calculations for both choices and show the results for the 3s\/2 orbit of 208 Pb in Fig. (9). 
The incident electron energy Ei = 500 MeV, the proton kinetic energy T p = 100 MeV, and 
the outgoing proton momentum is equal to the momentum transfer q. The solid line is 
the directly calculated R' LT and the dotted line is the apparent fourth structure function 
obtained when z is along q '(R). The dashed line is the directly calculated R' LT and the 
dash-dotted line is the apparent fourth structure function obtained by using q'(0) to define 
the z axis. 

When choosing the z axis along q'(0), the apparent structure function is out of phase with 
the direct structure function and the magnitude is suppressed at the first peak, but by using 
the z axis along q '(R) the structure functions are in phase and the magnitudes are quite 
close. Thus, changing the z axis to be along the direction of q '(R) permits the extraction of 
a fourth structure function R' LT . Furthermore, choosing this different z axis largely removes 
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the Coulomb distortion effects on the fourth structure function, at least around the first 
peak. From these results, we recommend that one can experimentally extract the fourth 
structure function by choosing the z axis along the new modified momentum transfer q '(R). 
Note that Op of the ejected proton in these plots is the polar angle measured from the 
differently chosen z axes. 

Using a polarized incident electron beam and detecting the knocked out proton out of 
the scattering plane, the fifth structure function R' LT , can be extracted by measuring the 
up-down asymmetry of the nuclear response. We also choose the incident electron energy 
Ei = 500 MeV, the proton kinectic energy Tp = 100 MeV, and the momentum transfer q 
equal to the proton momentum p for this case. 

We first look at the Coulomb distortion effect on the measurement of the fifth structure 
function from the 3 si/ 2 orbit of 208 P6 at fixed proton azimuthal angle <pp = 40° as shown 
Fig. (10). The direct fifth structure function again agrees with the apparent structure 
function in PWBA calculation as expected. The dotted line is the fifth structure function 
for the PWBA, and the solid line is the fifth structure function with the z axis along the 
asymptotic momentum q and the dashed line is for the case of the z axis along the modified 
momentum transfer q \R). We confirmed that the direct fifth structure function agrees 
with the apparent structure function for the approximate DWBA calculation as expected 
from our earlier observation that the approximate structure functions are symmetric under 
the transformation (f)p —» —(j)p. Thus when one calculates a u — a D all other terms cancel 
leaving only the R' LT i contribution. Note that unlike the case for R' LT , changing the z 
axis does not affect the shape and magnitude and does not reduce the Coulomb distortion 
effect significantly. Thus, one can experimentally extract a fifth structure function without 
redefining the z axis. Of course, Coulomb distortion clearly affects the magnitude ot R' LT , 
as compared to the plane wave result (by approximately 30% for 208 Pb). 

In Fig. (11), we show the fifth structure function as a function of <f>p at fixed polar proton 
angles, Op = 4°, for the 3 si/ 2 orbit of 208 Pb and the lsi /2 orbit of 16 0, and Op = 14° for the 
lpi /2 orbit of 16 0. These polar angles are the first peak position of the fifth structure function 
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for the S \/2 orbit of 208 Pb and the pi/ 2 orbit of 16 0 respectively for these kinematics. The 
dotted line is the PWBA result and the solid line and the dashed line are the approximate 
DWBA results obtained by choosing the z axis along the asymptotic momentum transfer q 
and along the momentum transfer q '(R). The approximate DWBA calculation for the pi / 2 
orbit of 16 0(e, e'p) in Fig. (11) reproduces the same shape as the full DWBA calculation [4] 
for the fifth structure function, but both differ in magnitude from the plane wave result. In 
this previous paper [4] where the full DWBA calculation was applied to a particular case, 
it was concluded that extracting the fifth structure function Rlt 1 at (j)p = 90° or averaging 
over (j)p largely removed the Coulomb distortion effects. Clearly that is not the case here. 
Therefore, it is not true in general that the Coulomb distortion effect can be removed at 
(j)p = 90°. Furthermore, we again note that choosing a different z axis does not help in 
removing the Coulomb effect for the fifth structure function unlike the case of the fourth 
structure function. We also confirmed that the fifth structure function extracted by using the 
incident electron helicity dependence or the up-down asymmetry agree exactly as expected 
from the —► —0p symmetry discussed earlier in the approximate Coulomb distorted form 
factors. 


V. CONCLUSIONS 

We have developed a plane-wave-like approximate solution to the Dirac equation in the 
presence of the static Coulomb held of a nucleus which agrees rather well with the exact 
partial wave solutions inside a sphere of approximately three times the nuclear radius. The 
limited spatial range of the approximation is not a serious restriction for electron induced 
nuclear processes since the bound state wavefunctions that enter any such process drop 
off exponentially outside the nuclear radius R. Using this approximate wavefunction, along 
with a few additional approximations, we also obtained an approximate DWBA potential 
valid for momentum transfers greater than about 350 MeV/c. This approximate potential 
has the same Dirac structure as the plane wave Mqllcr potential although it contains some 
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spatially dependent phase factors which destroy the azimuthal spatial symmetry about the 
momentum transfer direction q. The basic ingredients in our approximate potential are the 
static Coulomb potential of the target nucleus and the elastic scattering phase shifts for the 
incoming and outgoing electron energies in this potential. 

We have compared wavefunctions, potentials and cross sections for the (e, e'p) reaction 
on nuclei in the quasielastic region calculated with our approximation and previous approx¬ 
imations to the exact partial wave results. We find that the previous approximate results 
are in serious disagreement with the exact partial wave results. Our approximate results 
are in good agreement with the full DWBA results, but have a number of advantages over 
the DWBA results. The biggest advantage is that the plane-wave-like structure of the ap¬ 
proximate Coulomb distorted potential allows extraction of structure functions which are 
bi-linear products of transforms of the transition current components. However, unlike the 
PWBA analysis, these transforms are not just Fourier transforms, but contain additional 
spatial dependence on the kinematics resulting from the static Coulomb field of the target 
nucleus. 

In this paper, we investigated in some detail the extraction of the so-called fourth struc¬ 
ture function R' LT and fifth structure function R' LT , from the full cross section. We showed 
that this is possible, particularly for kinematics where the structure functions are large. 
However, for the case of the fourth structure function, R' LT) the z axis needs to be redefined 
to lie along the momentum transfer defined at the nuclear surface q \R) = p'(i?) — p 'f(R) 
to obtain agreement between the directly calculated structure function and the extracted 
structure function. For the fifth structure function R' LT ,, there is no need to redefine the z 
axis. 

The other major advantage of our approximate treatment of Coulomb distortion over the 
full DWBA partial wave calculation is that it is straightforward to apply it to higher energies. 
The full DWBA calculation at higher electron energies requires more and more partial waves 
with increasingly forbidding amounts of computer time needed. Using our approximate 
treatment, we do find it necessary to perform two additional numerical integrations over 
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angular coordinates 6 and 0 in the interaction matrix element as compared to a treatment 
that permits a multipole decomposition. However, these numerical integrations are not very 
time consuming and in a sense are just a replacement for summing over various intermediate 
angular momenta arising from angular momentum recoupling of the various partial wave 
expansions of the wavefunction and transition matrix elements. In our particular case, 
we have an analytic result for the electromagnetic potential so that the three dimensional 
integration is very fast as compared to a full partial wave analysis and the evaluation of 
thousands of radial matrix elements. 

In conclusion, our approximate treatment of Coulomb distortion for electron induced 
nuclear processes involving continuum nucleons works quite well and is particularly good 
for momentum transfers greater than 350 MeV/c. It permits the extraction of ’’structure 
functions” which should prove of great use in analysing (e,e'p) experiments at the Thomas 
Jefferson National Accelerator Facility and other laboratories. Of course, it is an approx¬ 
imation and if one has very high precision data it may be necessary to revert to the full 
DWBA calculation. 
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FIG. 1. Comparison of the exact scalar potential with approximate scalar potentials for the 
multipole component L=5, M=0. The initial energy Ej t = 400 MeV and the energy loss oj = 100 
MeV. The solid line is calculated with DWBA, the dotted line with EM A — k 2 and the dashed 
line is our approximate potential. 
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FIG. 3. Reduced cross sections for 20S Pb(e,e'p) from the 3 Si/ 2 shell with parallel kinematics. 
The kinematics are Ej = 412 MeV, and proton kinetic energy Tp = 100 MeV. The dotted line is the 
PWBA result and the dash-dotted line (which falls on top of it) is the approximate DWBA result 
with Z = 1. The dash-three-dotted line is the EM A — k 2 result, the solid line is the approximate 
DWBA result, the dashed line is the full DWBA result, and the diamonds are data from NIKHEF. 
The same spectroscopic factor of 71% is used in all curves. 
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FIG. 5. Reduced cross sections for 208 Pb(e,e'p) from the 3 Si/ 2 shell with perpendicular kine¬ 
matics. The kinematics are E t = 412 MeV, and proton kinetic energy Tp = 100 MeV. The dotted 
line is the PWBA result, the dash-three-dotted line is the EM A — n 2 result, the solid line is our 
approximate result, and the dashed line is the full DWBA result. 
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FIG. 7. Reduced cross sections for 208 P6(e, e'p ) from the 3^ / 2 shell for high missing momentum. 
The kinematics are E t = 487 MeV, momentum transfer q = 221 MeV/c, energy transfer lo = 110 
MeV, and proton kinetic energy Tp = 100 MeV. The solid line is the approximate DWBA result, 
the dotted line is the PWBA result, and the diamonds are data from NIKHEF. A previously 
determined spectroscopic factor of 71% was used for both curves. 
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FIG. 8. The fourth structure function for 208 P6(e, e'p) from the 3 s ^/2 shell as a function of 
missing momentum. The kinematics are Ej = 500 MeV, and proton kinetic energy Tp = 100 MeV. 
The dash-dotted line is the apparent structure function and dotted line is the directly calculated 
structure function for the PWBA (which fall on top of each other). The solid line is the apparent 
and the dashed line is the directly calculated structure function for the approximate DWBA result, 
while the dash-three-dotted line is the apparent structure function for the full DWBA result. 
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FIG. 9. The fourth structure function for 208 Pb(e,e'p) for the 3 Si/ 2 orbit as a function of the 
polar angle of the ejected proton. The solid line and the dotted line are the direct and the apparent 
structure function with the z axis along the momentum transfer q \R), and the dashed line and 
the dash-dotted line are the direct and the apparent structure function for the case of the z axis 
along the EMA momentum q'(0). 
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FIG. 10. The fifth structure function for 20S Pb(e,e , p) from the 3 Si/ 2 shell as a function of the 
polar angle of the ejected proton. The dotted line is the fifth structure function for the PWBA, 
while the solid and the dashed lines are approximate DWBA results with the z axis along the 
asymptotic momentum q and along the modified momentum transfer qd-R) respectively. 
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FIG. 11. The fifth structure function as a function of the azimuthal angle 4>p for 208 Pb(e,e'p) 
from the 2>S\/2 shell and for 16 0(e, e'p) from the lsi /2 shell and the lp \/2 shell. The dotted line 
is the fifth structure function for PWBA, the solid line is the fifth structure function with z axis 
along the asymptotic momentum q, and the dashed line is for the case of the z axis along the 
momentum transfer q '(R) . 
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